EUROPEAN ORGANISATION FOR NUCLEAR RESEARCH (CERN) 



CERN-EP/2000-052 
11 April 2000 



Measurement of W-pair production 
in e"^e~ collisions at 189 GeV 



The ALEPH Collaboration 



Abstract 

The production of W'''W~ pairs is analysed in a data sample collected by ALEPH at a 
mean centre-of-mass energy of 188.6 GeV, corresponding to an integrated luminosity of 
174.2 pb^^. Cross sections are given for different topologies of W decays into leptons or 
hadrons. Combining all final states and assuming Standard Model branching fractions, 
the total W-pair cross section is measured to be 15.71 ± 0.34(stat.) ± 0.18(syst.) pb. 
Using also the W-pair data samples collected by ALEPH at lower centre-of-mass 
energies, the decay branching fraction of the W boson into hadrons is measured to 
be B(W hadrons) = 66.97 ± 0.65(stat.) ± 0.32(syst.)%, allowing a determination of 
the CKM matrix element \Vcs\ = 0.951 ± 0.030(stat.) ± 0.015(syst.). 
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1 Introduction 



This letter presents results on W-pair production in e^e^collisions using data collected 
with the ALEPH detector at a centre-of-mass (CM) energy around 189 GeV, during the 
1998 data taking period. The WW events are identified in all possible W decay channels, 
thus allowing the determination of the W branching fractions and indirectly, the coupling 
of the W to cs pairs. 

The experimental conditions and data analysis follow those used in the cross section 
measurements at lower LEP2 energies. As they are already described in detail in |jl|], 
attention is focused here on changes to selection procedures other than a simple rescaling 
of cuts with the increased collision energy. 

A detailed description of the ALEPH detector can be found in Ref. and of 
its performance in Ref. 0. The luminosity is measured from small-angle Bhabha 
events, using lead-proportional wire sampling calorimeters [Q, with an accepted 
Bhabha cross section of approximately 4.25 nb f^. An integrated luminosity of 
174.20 ± 0.20 (stat.) ± 0.73 (syst.) pb""*^ was recorded at a mean CM energy of 
188.63 ± 0.04 GeV §. 

In this letter, the quoted signal cross sections are the CC03 cross sections 0, defined 
as the production of four-fermion final states through two resonating W bosons. Two 
processes contribute, z/g exchange in the t-channel and Z/7 exchange in the s-channel. 
The measured cross sections are corrected for the difference, denoted the "4f-CC03 
correction" [|I|, between the accepted cross sections for CC03 processes and all Standard 
Model four-fermion final states consistent with W-pair decays. 

The CC03 Standard Model cross section (cjww), calculated at \/s = 188.63 GeV with 
the program GENTLE g] is 16.65 pb (±2%). The KORALW version 1.21 Monte Carlo 
event generator is used to simulate the signal events with a normalised cross section in 
agreement with the GENTLE value. The JETSET []T0[ package is used for the hadronisation. 
Comparison samples, generated with EXCALIBUR [jri| and grc4f [0 for both CC03 and 
all four-fermion diagrams, are used for systematic error evaluation. Samples of events 
are also generated with different W masses, both for CC03 diagrams and for all WW-like 
four-fermion diagrams, with KORALW. 

The KORALZ [|13| Monte Carlo program is used to generate e"'"e~ — *• qq background 



events. Other backgrounds are generated with PYTHIA 5.7 flj] for ZZ, Zee and 



Weu processes, PH0T02 for 77 interactions, KORALZ for /i and r pair production 



and BHWIDE 161 and UNIBAB [171 for Bhabha events. 



2 Selection of W-pair candidates 
2.1 WW iuiu events 

The selection of fully leptonic W-pair decays follows exactly the two selections used for 
the cross section and branching fraction measurements at 183 GeV |T| . The two selections 
have similar overall efficiencies and background levels but differ in their sensitivities to the 
individual dilepton channels. The first is based on topological information and is sensitive 
to all channels. In the second, lepton identification is used to optimise the cuts according 
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to which final state is being considered. Events are accepted as WW candidates if they 
pass either of the two selections and are then classified into six di-lepton channels making 
use of electron and muon identification criteria. A jet or a single charged particle track 
is classified as a tau if no lepton is identified or the identified lepton has an energy less 
than 25 GeV. 

Beam related background, not simulated in Monte Carlo events, affects the efficiency 
of the cut which removes events depositing energy within 12° of the beam. Random 
trigger events were used to model these local energy deposits close to the beam. The 
inefficiency introduced is found to be 8.2 ± 0.5%. The CC03 efficiencies in the individual 
channels are given in Table |I| after correction for this beam related effect. The 
inclusive combination of the two selections has an overall efficiency of 64.2 ± 0.4% for the 
fully leptonic channels when combined assuming lepton universality. The total background 
amounts to 131 ± 7(stat.) ± 8.5(syst.) fb and is dominated by 77 — ^ U and non-WW-like 
ZZ — > events. In the data, the inclusive combination selects 220 events. 

All sources of systematic uncertainties are listed in Table 0. They are dominated by 
the uncertainties on the background cross sections, the uncertainty from the cut on energy 
detected close to the beam and by Monte Carlo statistics. 

A maximum likelihood fit is applied to determine the cross section for each fully 
leptonic decay channel using the efficiency matrices for signal and backgrounds given in 
Table ^ For all channels together the 4f-CC03 correction is — 10 ± 10 fb, where the 
uncertainty comes from Monte Carlo statistics. 

The results of the fit are 
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The systematic uncertainties are obtained by varying all input parameters in the fit 
according to their uncertainties. 

The total fully leptonic cross section is obtained with the same fit, assuming lepton 
universality: 

ct(WW ^ iuiu) = 1.78 ± 0.13(stat.) ± 0.02(syst.) pb, 
consistent with the sum of the individual channels. 

2.2 WW ii^qq events 

As for the lower energy measurements, three iuqq selection procedures are applied. One 
selection requires an identified electron or muon. The other two are designed to select 
rz/qq events, based on global variables or topological properties of the events. 

The selection of ei^qq and fiuqq events has been modified with respect to the previous 
analysis |jl| to take into account the greater initial boost of the W's. The preselection 
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Table 1: Summary of results of the different event selections on Monte Carlo and data events. Efficiencies 
are given in percent of CC03 processes; they assume Standard Model branching fractions in the overall 
values quoted in the right-hand column. The qqqq column refers to events with a NN output greater than 
0.3 where the backgrounds also include non-qqqq WW decays. The listed backgrounds do not include 
the 4f-CC03 corrections. 
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remains similar. It is based on the total charged particle energy and multiplicity and a cut 
on the longitudinal momentum and visible energy to reject Z7 events with an undetected 
photon. The selection of the lepton track relies on the fact that it is in general more 
energetic and isolated than the charged particles from the hadronic system. Thus the 
candidate lepton is chosen as the charged track that maximises ( 1 — cos Oej ) , where pg is 
the track momentum and 6ej the angle of the track to the closest of the jets clustered using 



the remaining reconstructed charged particles. For the jet clustering the DURHAM-P [|l^ 
algorithm with a Ucnt of 0.0003 is used. 

The same electron or muon identification criteria as for the fully leptonic channels are 
required for this lepton candidate track. However, no cut is applied on the lepton energy, 
so that Tuqq events where the r decays to a softer lepton (as r euu or r ^ ^uu) are 
also selected by this analysis. 

For electron candidates the lepton energy is corrected for possible bremsstrahlung 
photons detected in the electromagnetic calorimeter. The isolation of the lepton is defined 
as log(tan6'c/2) -|-log(tan 6*^/2) where 9c and 9p are, respectively, the angle of the lepton 
to the closest charged track, and the opening angle of the largest cone centred on the 
lepton direction which contains a total energy smaller than 5 GeV. 

For each event, probabilities that it comes from each of the three signal processes, 
euqq, /luqq or Tuqq, are determined (Fig. |I]) using Monte Carlo reference samples of 
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signal and backgrounds. These are evaluated from the identity, energy and isolation of 
the lepton plus the event total transverse momentum. An event is classified as euqq or 
/iz/qq if its corresponding probability is greater than 0.40; it is then not considered in the 
tau search. 

The selection of rz/qq events is based both on global event variables and a topological 
selection which attempts to identify the r jet. As the selections were described in detail in 
previous papers IQ , only changes other than a rescaling of the values of the energy-based 
cuts are described in the following: 

(i) in the global analysis the visible mass is required to be greater than 85 GeV/ and 
less than 155 GeV/ to take account of the boost of the W boson. The estimated 
energy of the "primary" neutrino must be smaller than 70 GeV; 

(ii) in the topological analysis, the energy of the most energetic quark jet must be less 
than 75 GeV and the mass of the hadronic system, i.e., the mass excluding the tau 
jet, is required to be less than 100 GeV/ c^. 

In addition, if an event with a well defined e or /i fails the e//i probability cut it is 
considered as a r candidate and kept if the ri/qq probability (Fig. [l|c) is greater than 0.4. 



Table 2: Systematic error summary (units in fb) 
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Figure 1: Probability distributions of preselected events for the (a) ei^qq, (b) /ii^qq and (c) ri/qq 
selections. The points are the data and the histograms the Monte Carlo predictions. The non-qq 
backgrounds include ZZ, Zee, Wei^ and r+r~ processes. 
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2.2.1 Results 



Table |l| gives the efficiencies for each selection. The inclusive efficiencies for the three 
semileptonic decay channels are 87.8 ± 0.4% for the electron channel, 91.6 ± 0.4% for the 
muon channel and 66.5 ± 0.5% for the tau channel, giving an 82 ± 0.3% average efficiency 
for WW ^ ii^qq, with a background of 314± 12(stat.) ±25(syst.) fb. The overall 4f-CC03 
correction is +8 it 40 fb. A total of 1066 events are selected in the data. 

The systematic uncertainties on the combined semileptonic cross section are 
summarised in Table 0. The largest contribution arises from the Monte Carlo statistical 
error on the 4f-CC03 correction. Uncertainties arising from the choice of the lepton 
isolation criterium and the probability cut are estimated from the change of efficiency 
following a bin-by-bin reweighting of the respective Monte Carlo one-dimensional 
distributions to the data. Background normalisation mainly affects the tau channel 
although there is also a contribution from residual Bhabha background in the electron 
channel. Decays of the Z to electrons and muons are used to evaluate lepton identification 
uncertainties whilst the contribution from beam related backgrounds is estimated by 
superimposing the energy deposits from random triggers on to the simulated events. 

To evaluate the individual cross sections a similar fit as for the fully leptonic events is 
used with the corresponding matrix of efficiencies and backgrounds. This yields 

(j(WW ^ ei/qq) = 2.41 ± 0.14(stat.) ± 0.07(syst.) pb, 
(t(WW ^ /ii/qq) = 2.39 ± 0.13(stat.) ± 0.06(syst.) pb, 
a(WW Tuqq) = 2.23 ± 0.17(stat.) ± 0.08(syst.) pb, 

where the systematic uncertainties are obtained by varying all input parameters of the fit 
according to their uncertainties. 

The total iuqq cross section is extracted by means of the same fit, under the 
assumption of lepton universality: 

a(WW ^ iuqq) = 7.07 ± 0.23(stat.) ± 0.12(syst.) pb, 
again consistent with a simple sum of the three channels. 



2.3 WW — > qqqq events 

The analysis of WW decays to four jets is updated from Ref. and consists of a simple 
preselection followed by a fit to the distribution of the output of a neural network (NN) 
with 14 input variables. 

In the preselection, the first step is to remove events with a large undetected initial 
state (ISR) photon from radiative returns to the Z by requiring that the modulus of the 
total longitudinal momentum of all objects is less than 1.5{Mms — Mz) where M^s is the 
observed visible mass. The particles are then forced into four jets using the DURHAM-PE 



algorithm ||18[ and the value of 1/34, where a four-jet event becomes a three-jet event, is 
required to be greater than 0.001. To reject qq events with a visible ISR photon none 
of the four jets can have more than 95% of electromagnetic energy in a one degree cone 
around any particle included in the jet. Four-fermion final states where one of the fermions 
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is a charged lepton are rejected by requiring that the maximum energy fraction of a single 
charged particle in a jet be smaller than 0.9. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

NN output 



Figure 2: Comparison of NN output distributions for data and Monte Carlo after the four-quark 
preselection. The points are the data and the histograms the Monte Carlo predictions. The line shows 
the fit result. 

At this point, 3438 events are selected in the data while 3593 ± 7 events are expected 
from Standard Model processes. This preselection has an efficiency of 98.4% for CC03 
events and a purity of 35.9%. 

The input variables for the NN are described in the Appendix and are related to the 
global event properties, the properties of jets, WW kinematics and the b-tag probabilities 
for the four jets. The NN is simplified with respect to that used for the analysis of the 
183 GeV data [0], with fewer input variables and a slightly improved performance. The 
output distributions of the NN for the data compared with the signal and backgrounds 
predicted by the Monte Carlo are given in Fig. |^. In Table |l| results are given for a cut 
at 0.3 on the NN value. 

The cross section is extracted by means of a binned maximum likelihood fit to the 
full NN output distribution. In the fit only the normalisation of the Monte Carlo signal 
is allowed to vary; all backgrounds, including WW £z/qq, are kept fixed both in shape 
and normalisation. Systematic studies at the Z peak when two jet qq events are forced 
into four jets show small discrepancies between data and Monte Carlo. This affects the 
qq background shape and is quantified using the 14 NN variables to yield a correction of 
+50 ± 50 fb to the value from the fit. The corrected fit result is: 

a(WW ^ qqqq) = 6.89 ± 0.23(stat.) ± 0.13(syst.) pb. 
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If only events with a NN output greater than 0.3 are kept a value of 6.89 ± 0.23 pb is 
found showing no significant bias due to the qq background. 

The sources contributing to the systematic uncertainty are summarised in Table |^. 
The largest contribution is the uncertainty in the qq background normalisation where 
a 5% variation is assumed. The contributions from uncertainties in the values of jet 
variables used in the event preselection and the NN input are assessed by adjusting their 
directions and energies according to residual differences observed between Z calibration 
data and Monte Carlo fO. Possible miscalibrations of the calorimeters are also taken into 



account. The beam background is studied in the same way as described in section |2.2.1 
The Monte Carlo statistics error is dominated by the 4f-CC03 correction. 

Uncertainties in the WW generator are evaluated by comparing samples of fully 
simulated Monte Carlo events generated with KORALW and EXCALIBUR. To establish a 
WW fragmentation uncertainty, the HERWIG |2^ generator was tuned at the Z peak both 
for all flavours and non-b quark flavors. Then, the same samples of signal events generated 
with KORALW are fragmented with both JETSET and the appropriately tuned HERWIG. For 
the qq background uncertainty, a sample of events generated with KORALZ using JETSET is 
compared with a separate sample of pure HERWIG events to assess the effect of the choice of 
generator. Colour reconnection effects are estimated using the SKI model in JETSET with 
a reconnection probability of 0.3 and the effect of Bose-Einstein correlations is estimated 
according to the scheme, denoted BE3, as proposed for the LUBOEI |2ll] implementation 
in JETSET. The procedures followed are the same as those used in Ref. |19[ . 

Several cross checks have been performed on the fit result to search for possible biases 
arising in the selection and full simulation of the WW events. The previous version 
of the analysis with a different preselection and NN [|l| gives 6.81 ± 0.23 pb. Another 
estimate with a preselection based mostly on charged tracks and a six variable NN using 
only charged tracks gives 6.83 ± 0.25 pb. Also, a selection based purely on calorimeter 
measurements with six input variables to a linear discriminant gives 6.91 ± 0.26 pb. A 
variety of linear discriminant analyses using from 4 to 14 variables give results which vary 
from 6.65 to 6.77 pb. All these checks give results which are consistently lower than the 
GENTLE prediction. 



3 Total cross section 

The total cross section is obtained from a fit to all channels described above assuming the 
Standard Model branching fractions, the only unknown being the total cross section. The 
fit uses the matrices of efficiencies and backgrounds for the various analyses and yields 

dww = 15.71 ± 0.34(stat.) ± 0.18(syst.) pb. 

Assuming no additional unexpected decay mode, the result is not significantly different if 
the branching fractions of the Standard Model decay modes are unconstrained. 

The measurement is 5.5% lower than the GENTLE prediction. However, new calculations 
including full 0{a) electroweak corrections, calculable in the double pole approximation 



(DPA) [^] have recently appeared. Two Monte Carlo programs, YFSWW3 p3| and 



RacoonWW p^], are being developed. First numerical calculations find cross sections. 
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respectively, 1.9% [g3[ and 2.4% lower than GENTLE at 189 GeV. The predictions 
from RacoonWW also include soft-photon exponentiation and leading log corrections for 



initial state radiation beyond 0{a) in addition to the calculations described in Ref. |24 . 
The uncertainty in the two new models is expected to be of the order of 0.5% 

Fig. 1^ shows the total cross section measured as a function of the CM energy. The 
predictions of the two more complete YFSWW3 and RacoonWW calculations are also shown 
and are in better agreement with the experimental results than GENTLE. At 189 GeV, the 
measurement is 3.8% (1.5 standard deviations) below the YFSWW prediction and 3.2% 
(1.3 standard deviations) below the RacoonWW prediction. Taking into account also 
the cross section values measured by ALEPH |l| at 172 and 183 GeV, the data are 
4.6 ±2.0%, 2.8 ±2.0%, and 2.3 ±2.0% below the GENTLE, YFSWW and RacoonWW predictions 
respectively. These values use the signal efficiencies determined with KORALW and the 
quoted systematic uncertainty takes no account of any efficiency difference which may 
arise from the new calculations. 
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Figure 3: Measurements of the W-pair production cross section at four CM energies, compared with 
the Standard Model predictions from GENTLE, YFSWW3 and RacoonWW, for the LEP average value of the 
W mass [Eel. 
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4 Branching fractions and Vcs 

The same fit as for tlie total cross section is performed combining the data samples 
collected at 161, 172, 183 and 189 GeV CM energies. 

Without assuming lepton coupling universality, the seven unknowns are the three 
individual leptonic branching fractions and the four total cross sections at 161, 172, 183 
and 189 GeV. The hadronic branching fraction is set to 1 — Be — — B^-. The fitted 
leptonic branching fractions are 

B(W^ei/) = 11.35 ± 0.46(stat.) ±0.17(syst.)%, 
B(W-^/ii^) = 11.10 ± 0.44(stat.) ±0.16(syst.)%, 
B(W^rz/) = 10.51 ±0.55(stat.) ±0.22(syst.)%, 

and are consistent with lepton universality and the Standard Model expectations. Due to 
cross-contaminations in the identification of W decays to rz/, eu or fiu, the measured 
B(W Tu) is 26% anticorrelated with B(W eu) and 24% anticorrelated with 
B(W fxu). The B(W ^ eu) is 4.6% anticorrelated with B(W ^ fxu). 

If lepton universality is assumed a fit for B(W qq) and the total cross sections at 
the four energies yields 

B(W ^ qq) = 66.97 ± 0.65(stat.) ± 0.32(syst.)%. 

This result can be expressed in terms of the individual couplings of the W to quark- 
antiquark pairs: 

B(W^qq) 2 , IT/ 12 , IT/ 12 , IT/ 12 , IT/ 12 , |t/ |2^/1 , ^, /_2 



B(W-^qq) " (lKdr+|Kdr + |Ksr+|Ksr + |Kbr + |Kbr)(l + as(m2v)/vr). 



:ii 



The least well known of these is |V^s|- Using the world average value of as(m|) evolved to 
as(^w) — 0.121 ± 0.002, and the squared sum of the other measured CKM matrix 



W ; 



elements which is 1.05 ± 0.01, the measured hadronic branching fraction is 

\Vcs\ = 0.951 ± 0.030(stat.) ± 0.015(syst.). 



5 Conclusions 

The W-pair production cross section at -y/i = 188.63 GeV has been measured in all decay 
channels from an integrated luminosity of 174.20 pb'^. The total cross section is found 
to be 

aww = 15.71 ± 0.34(stat.) ± 0.18(syst.) pb. 

This result is 5.5% (1.9 standard deviations) lower than the GENTLE prediction but 
in better agreement with more recent calculations. It agrees well with the recent 



measurement by the DELPHI Collaboration p8| at the same CM energy. 

After inclusion of the data taken at CM energies of 161, 172 and 183 GeV, the hadronic 
decay branching fraction is found to be 66.97 ± 0.65(stat.) ± 0.32(syst.)% which is used 
to determine the CKM matrix element \Vcs\ equal to 0.951 ± 0.030(stat.) ± 0.015(syst.). 
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Appendix: Neural network input variables for the qqqq 
selection 

The neural network hadronic event selection uses 14 variables. These are based on global 
event properties, heavy quark flavour tagging, jet properties and WW kinematics and are 
listed below. The four jets are numbered in order of decreasing energy. 

Global event properties 

• Thrust 

• Sphericity 

• Missing energy 

• Sum of the four smallest intcrjct angles 
Heavy flavour tagging 

• Probability of an event being a light quark (uds) event based upon impact parameter 
significance of charged particles in the event 

Jet properties 

• Maximum electromagnetic energy fraction of a jet in any one degree cone 

• Maximum summed charged particle energy fraction of a jet 

• Minimum number of charged particles in a jet 
WW kinematics 

• Angle between Jet2 and Jet3 

The following jet related variables are determined from kinematically fitted jet 
momenta. 

• Energy of Jetl 

• Energy of Jet3 

• Energy of Jet4 

• Smallest jet mass 

• Second smallest jet mass 
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